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Abstract 

Understanding the origin of high-mass stars is central to modern astrophysics. We shed light 
on this problem with simulations using a novel, adaptive-mesh, ray-tracing algorithm. These 
simulations consistently follow the gravitational collapse of a massive molecular cloud core, 
the subsequent build-up and fragmentation of the accretion disk surrounding the nascent 
star, and, for the first time, the interaction between its intense UV radiation field and the 
infalling material. We show that ionization feedback can neither stop protostellar mass 
growth nor suppress fragmentation. We discuss the effects of feedback by ionizing and 
non-ionizing radiation on the evolution of the stellar cluster. The accretion is not limited 
by radiative feedback but by the formation of low-mass companions in a process we call 
"fragmentation-induced starvation" . This behavior consistently reproduces the observed 
relation between the most massive star and the total mass of stars in a cluster. We show 
that magnetic fields reduce the star formation rate and lead to the formation of more massive 
stars. 



1 Introduction 



High-mass stars form in denser and more massive cloud cores (Motte et al. , 2008) than their 



low-mass counterparts (Myers et al. , 1986). High densities also result in local gravitational 



instabilities in the accretion flow, resulting in the formation of multiple additional stars 



(Klessen & Burkert, 2000 Kratter & Matzner 



observed to have companions (Ho &; Haschick 



2006 ) . Young massive stars are almost always 
81), and the number of their companions 
significantly exceeds those of low- mass stars (Zinnecker &: Yorke, 2007). Such companions 
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influence subsequent accretion onto the initial star (Krumholz et al. , 2009). Observations 



show an upper mass limit of about 100 Mq. It remains unclear whether limits on internal 
stability or termination of accretion by stellar feedback determines the value of the upper 
mass limit (Zinnecker & Yorke, 2007). 

The most significant differences between massive star formation and low-mass star 
formation seem to be the clustered nature of star formation in dense accretion flows and the 
ionization of these flows. We present the first three-dimensional simulations of the collapse 
of a molecular cloud to form a cluster of massive stars that include ionization feedback by 
Peters et al. ( 2010a[c|b 2011), allowing us to study these effects simultaneously. 



2 Numerical Method and Initial Conditions 



We present the first three-dimensional, radiation-(magneto-)hydrodynamical simulations of 
massive star formation, taking into account heating by both ionizing and non-ionizing radi- 
ation, using the adaptive- mesh code FLASH (Fryxell et al. , 2000). We propagate the radia- 
tion on the adaptive mesh with our extended version of the hybrid characteristics raytracing 



method (Rijkhorst et al. , 2006 Peters et al. , 2010a). We use sink particles (Federrath et al. 



2010) to model young stars. Sink particles are inserted when the Jeans length of collapsing 



gas can no longer be resolved on the adaptive mesh. They continue to accrete any high- 
density gas lying within their accretion radius. We use the sink particle mass and accretion 
rate to determine the radiation feedback with a prestellar model (Peters et al. , 2010a). A 
detailed description of our numerical method and the underlying assumptions can be found 
in 



Peters et al. (2010a). 



We start our simulations with a 1000 Mq molecular cloud having a constant density 
core with p = 1.27 x 10"'^'^ gcm~^ within a radius of r = 0.5 pc, surrounded by an r~^/^ 
density fall-off out to r = 1.6 pc. The cloud rotates as a solid body with an angular velocity 
Lo = 1.5 X 10~^^s~^. The initial temperature is T = 30 K. The highest resolution cells on 
our adaptive mesh have a size of 98 AU. Sink particles are inserted at a cut-off density of 
/Jcrit = 7 X 10~^^ gcm~^ and have an accretion radius of rsink = 590 AU. 

We compare the results of four different simulations. In the first simulation (Run A), 
a dynamical temperature fioor is introduced to suppress secondary fragmentation. Only one 
sink particle (representing a massive protostar) is allowed to form. In the second simulation 
(Run B), secondary fragmentation is allowed, and many sink particles form, representing a 
group of stars, each contributing to the radiative feedback. The third simulation (Run D) 
is a control run in which secondary fragmentation is still allowed, but no radiation feedback 
is included. The fourth simulation (Run E) is a magnetized version of Run B, the full 
stellar group simulation with radiation feedback from all stars. Run E includes an initially 
homogeneous magnetic field along the rotation axis of the cloud with a magnitude of 10 //G, 
corresponding to a mass-to-fiux ratio (M/$) = 14(M/$)cr in the central core at the beginning 
of the simulation. See Peters et al. (2010b) and Peters et al. (2011) for a thorough discussion 
of these initial conditions. 
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Figure 1: Accretion history of the single (Run A) and multiple (Run B) sink simulations. 
Run A was stopped at 72 M©, while no sink particle in Run B exceeds 25 Mq over the 
simulation runtime. The left hand plot shows the sink particle masses for Run A (green), 
the individual sink masses of Run B (blue) as well as the total mass in sink particles in Run 
B (red). The right hand plot shows the accretion rates of the sink particle in Run A (green) 
and the sink particle which forms first in Run B (blue), which also turns out to end up as the 
most massive at the end of the simulation. While the accretion rate in Run A never drops 
below 10"^ Mq yr ^, accretion onto the most massive sink can drop significantly below this 
value and even be stopped totally in Run B. 



3 Fragmentation-Induced Starvation 

We first examine the accretion histories of our different models. Figure [T] shows that in 
Run A, with only one sink particle allowed to form, nothing halts accretion onto the central 
protostar. It continues to grow at an average rate of M ~ 5.9 x 10~^ Mq yr~^ until we stop 
the calculation when the star has reached 72 Mq. The increasingly massive star ionizes the 
surrounding gas, raising it to high pressure. This gas breaks out above and below the disk 
plane, but it cannot halt mass growth through the disl^ mid-plane. 

In contrast, in Run B, where multiple sink particles are allowed to form, two subsequent 
sink particles form and begin accreting soon after the first one, and many more follow within 
the next 10^ yr (see Figure [l]). By that time the first sink has accreted 8Mq. Within another 
3 X 10^ yr seven further fragments have formed, with masses ranging from 0.3 Mq to 4.4 Mq 
while the first three sink particles have masses between 10 Mq and 20 Mq, all within a radius 
of 0.1 pc from the most massive object. Accretion by the secondary sinks terminates the 
mass growth of the central massive sinks. Material that moves inwards through the disk 



driven by gravitational torques accretes preferentially onto the sinks at larger radii (Bate 



2002). Eventually, hardly any gas makes it all the way to the center to fall onto the most 
massive objects. This fragmentation-induced starvation prevents any star from reaching a 
mass > 25 Mq in this case. 



^We will hereafter refer to the flattened, dense, accretion flow that forms in the midplane of our rotating 
core as a disk. However, it is not necessarily a true Keplerian, viscous, accretion disk, which probably only 
forms within the central few hundred astronomical units, unresolved by our models. 
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Figure 2: (Left) Total accretion history of all sink particles combined forming in Runs A, B, 
and D. While the heating by non-ionizing radiation does not affect the total star formation 
rate, the ionizing radiation appreciably reduces the total rate at which gas converts into stars 
once the most massive object has stopped accreting and its H ii region can freely expand. The 
slope of the total accretion history in Run B goes down because the massive stars (dashed 
line) accrete at a decreased rate, while the low-mass stars (dotted line) keep accreting at the 
same rate. (Right) The total mass in sink particles Mginks ^is a function of the most massive 
star in the cluster Mmax- We plot the curves for Run A, Run B and Run D as well as the 
fit from the competitive accretion simulations by Bonnell et al. (2004). The turn-off away 
from the scaling relation (indicated by arrows) is shifted towards higher masses by radiative 
feedback. 



Figure [T] also reveals that the most massive stars in the cluster are those which form 
first and then keep accreting at a high rate. The most massive stars in Run B are the same 
for all times, but their mass keeps increasing during the formation of the stellar cluster. 

The accretion behavior in Run B contrasts sharply with competitive accretion models 
(Bonnell et al. , 2001, 2004), which have no mechanism to turn off accretion onto the most 
massive stars. In these models, material falls all the way down to the massive stars sitting in 
the center of the gravitational potential which thereby take away the gas from the surround- 
ing low-mass stars. In our fragmentation-induced starvation scenario, exactly the opposite 
happens. Figure [2] illustrates that, although the accretion rates of the most massive stars 
(M > 10 Mq) steadily decrease, the low-mass stars (M < 10 Mq), keep accreting at the same 
rate. 

Competitive accretion models show a correlation between the mass of the most massive 
star Mjnax and t he total cluster mass Mginks during the whole cluster evolution that is roughly 
Afmax oc Mgjj^j^g (Bonnell et al. , 2004). This correlation has been argued to represent a way to 
observationally confirm competitive accretion (Krumholz & Bonnell, 2009) and is in fact in 



good agreement with observations (Weidner & Kroupa, 2006; Weidner et al. , 2010). However, 
we find that our simulations also reproduce the observed relation between Mmax and Mginks- 
Figure [2] shows Mginks as function of Mmax for Run A, Run B and Run D, and the 



relation Mn 



0.39M 



2/3 



sinks' 



which was found by 



Bonnell et al. 



(2004) as a fit to their 



simulation data. Over the whole cluster evolution, the curve for Run D lies above this fit. 
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Figure 3: Individual accretion histories for Run B and Run D. The figure shows the stellar 
masses as function of time for all sink particles that form in Run B (left) and Run D (right). 
Because the Jeans mass is lower without radiative feedback (Run D), many more sink particles 
form in Run D than in Run B, and the mass of the most massive stars is also lower. 



while the curve for Run B always lies below it. The fit agrees with our simulation data as 
well as it does to that of Bonnell et al. (2004). This indicates that the scaling is not unique 
to competitive accretion, but can also be found with the fragmentation-induced starvation 
scenario and hence cannot be used as an observational confirmation of competitive accretion 
models. 



4 Effects of Radiative Feedback and Magnetic Fields 

Figure |3] shows the individual accretion histories of each of the sink particles in Run B and 
Run D. Radiative heating cannot prevent disk fragmentation but raises the local Jeans mass. 
Hence, many fewer stars form in Run B than in Run D, and the mass of the most massive 
stars in Run B is higher. It is also evident from the figure that star formation is much more 
intermittent in the case with radiation feedback (Run B). The reason for this behavior is 
that the star formation process is controlled by the local Jeans mass, which depends to a 
large degree on how the filaments in the disk shield the radiation. Shielding can lower the 
Jeans mass temporarily and thereby allow gravitational collapse that would not have occurred 
otherwise. 

Since the accretion heating raises the Jeans mass and length in Run B, the total number 
of sink particles is higher in Run D than in Run B, and the stars in Run D generally reach 
a lower mass than in Run B. These two effects cancel out to lead to the same overall star 
formation rate for some time (see Figure [2]). At one point in the evolution, however, also 
the total accretion rate of Run B drops below that of Run D. At time t ~ 0.68 Myr the 
accretion flow around the most massive star has attenuated below the value required to trap 
the H II region. It is able to break out and affect a significant fraction of the disk area. 
A comparison with the mass growth of Run D clearly shows that there is still enough gas 
available to continue constant cluster growth for another 50 kyr or longer, but the gas can no 
longer collapse in Run B. Instead, it is swept up in a shell surrounding the expanding H II 
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Figure 4: (Left) Total accretion history of Run A, Run B, Run D and Run E. While the 
total accretion rates of Run B and Run D agree until ionization feedback stops massive star 
formation in Run B, the total accretion rate in Run E already starts to decline after 20kyr. 
The magnetic field in Run E additionally supports the gas against collapse, reducing the total 
accretion rate. (Right) Protostellar masses of the first sink particles in Run A, Run B, Run D 
and Run E. Radiative heating (Run B) and presence of magnetic fields (Run E) increase the 
final masses of the massive stars. The largest part of the additional mass accretion in Run E 
compared to Run B is due to the stronger initial accretion phase. 



region. 

It is also notable that the expanding ionization front around the most massive stars 
does not trigger any secondary star formation, which suggests that triggered star formation 
(Elmegreen &; Lada, 1977) may not be as efficient as expected, at least on the scales considered 
here. 

The total accretion histories of all sink particles combined in each of the four simulations 
are contrasted in Figure [4} For the first 20kyr, the total accretion rate of Run B, Run D and 
Run E is nearly identical. The accretion rate in these multiple sink simulations is generally 
higher than in the single sink Run A since the large group of sinks accretes from a large 
volume, without needing to rely on outward angular momentum transport to deliver material 
to the direct feeding zone of the central sink. In the magnetized Run E, no secondary sink 
particles form during the first 20kyr, so that the increased accretion rate in this phase is due 
to the additional angular momentum transport performed by the magnetic field. After the 
initial 20kyr, the total accretion rate in Run E falls below Run B and the control Run D, 
deviating increasingly with time, but always staying above the accretion rate in Run A. The 
accretion histories of Run B and Run D only start to separate at relatively late time when 
the ionizing radiation begins to terminate accretion onto the most massive sinks in Run B 
(see the discussion in Peters et al. 2010b), but the magnetic field additionally reduces the 
rate at which gas collapses in Run E. Thus, the total accretion rate with magnetic fields is 
lower than without magnetic fields. 

The central star in Run E can grow to larger masses because in the initial phase 
fragmentation is delayed by magnetic support. This allows the central object to maintain a 
high accretion rate for a longer time. Nevertheless, the accretion rate of the massive star in 
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Run E drops significantly when secondary sink particles form, since they form in a dense ring 
around the central massive sink and to some extent starve it of material, even though they 
never cut off accretion entirely. 

The stronger initial accretion phase in Run E compared to Run B yields the main 
contribution to the larger final mass of the massive sink. The magnetic field very efficiently 
redistributes angular momentum, resulting in an increased radial mass flux through the high- 
density equatorial plane. Consequently, the initial accretion rate in the magnetized simulation 
(Run E) even lies considerably above the non- magnetized single sink calculation (Run A). 



5 Conclusions 

We have reviewed some of the results of our recent radiation-(magneto-)hydrodynamical 
simulations of massive star formation (Peters et al. , 2010a|c|b 2011[ ), which for the first time 



simultaneously include the effect of heating by both ionizing and non-ionizing radiation. We 
find that ionization feedback is unable to stop protostellar mass growth. Instead, the mass 
of the most massive stars is limited by the formation of lower-mass companions in their 
gravitationally unstable accretion flow in a process we call fragmentation-induced starvation. 
Our numerical model reproduces the observed relation between the most massive star in a 
cluster and the total cluster mass. We find that heating by non-ionizing radiation decreases 
the total number of stars formed and increases the mass of the most massive stars, but does 
not change the total star formation rate. The star formation rate is reduced by ionizing 
radiation once the H ii regions can steadily expand. Even initially very weak magnetic 
fields can markedly reduce the star formation rate by providing additional support against 
gravitational collapse and lead to the formation of more massive stars via magnetic braking 
of the accretion flow. 
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